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Abstract 

We investigate a class of operators with non-vanishing averages in a D- 
dimensional matrix model recently proposed by Kazakov and Migdal. Among 
the operators considered are "filled Wilson loops" which are the most reason- 
able counterparts of Wilson loops in the conventional Wilson formulation of 
lattice QCD. The averages of interest are represented as partition functions of 
certain 2-dimensional statistical systems with nearest neighbor interactions. 
The "string tension" a', which is the exponent in the area law for the "filled 
Wilson loop" is equal to the free energy density of the corresponding statisti- 
cal system. The continuum limit of the Kazakov-Migdal model corresponds 
to the critical point of this statistical system. We argue that in the large N 
limit this critical point occurs at zero temperature. In this case we express 
a' in terms of the distribution density of eigenvalues of the matrix-valued 
master field. We show that the properties of the continuum limit and the 
description of how this limit is approached is very unusual and differs drasti- 
cally from what occurs in both the Wilson theory {S oc (Tr H U + c.c.)) and 
in the "adjoint" theory {S oc iTrHf^P)- Instead, the continuum limit of the 
model appears to be intriguingly similar to a c > 1 string theory. 

^This work is supported in part by the Natural Sciences and Engineering Research 
Council of Canada 

^Permanent Address: ITEP, 117295, Moscow, Russia. 



1 Introduction 



In a recent paper, |jT[| V.Kazakov and A.Migdal proposed a new lattice gauge 
model where the Yang-Mills interactions are induced by minimal coupling 
to a scalar field which transforms in the adjoint representation of the gauge 
group. Their theory (to be referred below as the Kazakov-Migdal model 
(KMM)) resembles a D-dimensional matrix model with partition function of 
the form ^ 

Zkmm= j[dU]d^eM-T.^^V{^x)) + 



ir^{x)U{x,y)^{y)U\x,y)} (1) 

<x,y> 

Here $ and U are N x N Hermitean and unitary matrices defined on the sites 
and the links of a D-dimensional rectangular lattice respectively. Although 
the potential V{^) is quite general, for the purpose of most of our work we 
can use the quadratic potential: 

V{^) = m^$2 (2) 



Throughout this paper we use square brackets to denote non-trivial measures of 
integration. Thus, for example, dX = W.i j dXij, but the Haar measure [dU] is not simply 
Yii i dUij, there is also a non-trivial factor, see (35). In fact [dU] will always denote the 



Haar measure over unitary matrices. We shall denote normalized averages in the theory 
(0) by <C while unnormalized integrals over unitary matrices with Haar measure [dU] 
will be denoted by < > . Unitary matrices U are taken to be elements of the group 
SU{N). If we were to allow U to be in U{N) then the Zjy-gauge invariance of the KMM 
to be discussed throughout the text would extended to a J7(l)-gauge invariance. This does 
not influence our results, since the model (|l|) is insensitive to the U{1) factor in U. 

Hermitian matrices such as $ and are not assumed to be traceless. However if the 
potential is quadratic the trace part of <i>(x) {^{x) — 4>ti-{x) ■ I + traceless part ) is easily 
separated from the traceless degrees of freedom in giving rise to an extra factor 



n 



The critical value of the bare mass for this abelian subsector of the KMM is obviously 
m^j.;j — D. This is the the same value of TO^it as is obtained for the non- abelian compo- 
nents of <i> in the quasi-classical approximation. 
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The effects of non-quadratic terms in the potential can be encoded in quan- 
tities such as the distribution function, of the eigenvalues of 

For D = 4 this model is a lattice approximation to a Yang-Mills theory 
alternate to the conventional Wilson lattice gauge theory. Recall that the 
partition function of the Wilson theory is given by 

/[rft/]exp(^X:^W^(°) + c-c.), (3) 

where IV(n) is the elementary Wilson loop, i.e. the trace of the product of 
link operators, around an elementary plaquette which we denote by □. 

The main advantage of the KMM is the existence of two reformulations of 
it, which, while being equivalent, look very different and emphasize different 
features of the theory. 

One reformulation arises from integration over the $ variables in (|T]) and 
can be written as 

Zkmm = J Mf/]e-^^//[^] cx / [dU] exp Q ^(^^2^(0 1^(^)1') , (4) 
where the sum in the action is over all oriented contours F on the lattice and 
l^(F)=tr( n U{x,y)]- (5) 

with 

W{-T) = W(T)- U{y, x) = W{x, y). (6) 

In contrast to the Wilson theory, the effective action depends on the modulus 
of VF(F), rather than on its real part. It also depends on contours of all sizes 
rather than on just the elementary plaquette. 

Another reformulation of the KMM arises by integrating first over U 
rather than over $. To do this, it is convenient also to integrate out the angu- 
lar components of $ leaving an integral over its eigenvalues 0i(x),...,0Ar(x) . 
The partition function is then given by 

Zkmm oc / WWdU^) e-^(*'(^» IT I[mA{y)] (7) 

X i <x,y> 
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where 



i<j 



and 



I[(f),iP] = / [dU] exp (trmW^) 



(8) 



The first representation (^) has been used to give naive arguments that the 
KMM (P is related to QCD (or, more precisely to quantum (7/Modynamics) 
if D=4. This idea, which we shall further refer to as the "naive continuum 
limit", is that the action in (^ clearly has an absolute minimum at U = I. 
Assuming that in the continuum limit the dominant contribution to the par- 
tition function comes from the vicinity of this minimum the effective action 
can be computed by writing U = exp{iaA) (where a is the lattice spacing). 
The result is derived in Ref. [||]: 



We shall argue in this paper (see especially Section 6) that this argument 
is too naive and that the approach to the continuum limit is much more 
complicated. 

The second representation (0-§) can be used for a much more reliable 
investigation of the KMM. It may even be possible to solve the KMM ex- 
actly for any N (although most of the progress has been made for large N). 
The reason for these hopes is that the integrals such as (P) have a pure al- 
gebraic interpretation so that standard methods such as Fourier analysis on 
co-adjoint orbits or the Duistermaat-Heckmann (DH) integration formalism 
can be applied to them. 

As argued in the transformation from (|^) to (0) is highly non-trivial 
from a conceptual point of view and is more or less equivalent to a summation 
of all planar diagrams in the Yang-Mills theory. 





with 




(10) 
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In this paper we report some preliminary considerations of correlation 
functions in the KMM which should correspond to the physical observables 
of "induced QCD". We shall concentrate our analysis on finite values of N. 
We shall then present arguments that the features of the model for finite N 
are sufficiently rich and interesting that they may have implications for the 
case of infinite N. 

We shall see that the KMM differs from the conventional Wilson theory 

in many respects but in particular in the properties of the continuum 
limit. We show in Section 6 that the probability distribution for the func- 
tional integral near U = I is much less sharp in the KMM than it is in the 
Wilson theory. Moreover the continuum limit in the KMM appears sensitive 
to the quantum measure in functional integral, not only to the form of the 
action. This has a crucial effect on the normalization of observable operators 
and probably on the other features of continuum limit as well. 

Another purpose of this paper is to develop a formalism for evaluation 
of the expectation values of operator observables in the KMM. We shall 
see that the most interesting observables in the KMM are associated with 
2-dimensional sublattices S of the original D-dimensional lattice. Vacuum 
expectation values of these quantities can be represented as partition func- 
tions of 2d statistical systems defined on S. We shall see that the Boltz- 
mann factors on the links of this system are fiuctuating variables so that 
generically they behave in a manner similar to spin glasses. In the zero- 
temperature limit, however, they turn into ordered magnetic-type systems 
similar to A^-state Potts models. We shall argue that in the limit of large N 
the zero-temperature limit of the 2d system can be identified with the con- 
tinuum limit of the KMM. Furthermore we can show using these ideas that 
the area law for filled Wilson Loops holds in this limit. The string tension 
a' can be expressed in terms of the density p(0) of eigenvalues of the master 
field (See Sections 5 and 7 for details.) 

Note that there is an "adjoint" model "in between" the KMM and Wilson 
theory. Consider a model with a partition function 

Za= [[dU]expY.\\tTWiO)f (11) 

Many of our comments in Sections 2 and 4 about the choice and the properties 
of observables in the KMM are a result of gauge invariance of the theory. 
They are thus equally true for Eq. (|1T|). The unusual behaviour of the 
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continuum limit which will be discussed in Sections 6 and 7 is peculiar to 
the KMM and is drastically different from that in the theories (|3| JTT|) . The 
reason for this difference is that in (^) there is a sum over all contours (not 
only over elementary plaquettes). In any case one must keep firmly in mind 
that the KMM is significantly different from both the Wilson theory from 
the model of Eq. (|lD. 

2 Correlators and Zjv- invar iance 

Any lattice gauge theory possesses a gauge invariance with respect to trans- 
formations specified by P{x) G U {N) at every site x and which act on the $ 
and ?7-fields by 

^{x) ^ P{x)^{x)P\x), (12) 

U{x,y)^P{x)U{x,y)P^{y). (13) 

The gauge invariant observables in lattice gluodynamics are made from traces 
of products of [/-matrices along closed oriented contours F on the lattice, i.e. 
conventional Wilson loops (|^,|^). 

The first important property of the KMM (1.1), which makes it different 
from Wilson's theory, is that the space of gauge-invariant operators must be 
further reduced P]. The reason for this is the presence of an additional Z{N) 
gauge symmetry: the freedom of independent discrete transformations 

U (x, y) uj{x, y)U{x, y); (14) 
^{x,y)^ = 1, 

so that LV is any A^-th root of unity. This additional gauge invariance excludes 
the conventional Wilson loop Vr(r) from the set of gauge invariant observ- 
ables (so that <^ W(T) 3>= 0). In fact this symmetry excludes any product 
of Wilson loops as an observable in the theory unless every link involved is 
passed an even number of times, equally in one and in the other direction. 
To be precise, the difference between the number of f/'s and f/^'s at every 
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link should be divisible by A^. The set of gauge invariant observables in the 
KMM thus consists of the products 

W{T,, r,) = j^W{T,)W{T2)...W{Tk), (15) 

provided that 

Ti + ... + rfc = [modA^] (16) 

(The contours are, of course, taken to be oriented). 

The simplest examples of allowed operators are the "Adjoint Wilson 
Loop" (Fig.l) 

W{T,-T) = ^\W{T)\^ (17) 

the "Baryon Loop" 

Wbaryon(r) = j^W (NT) (18) 

(in which contour wraps A^ times around F) and the "Filled Wilson Loop" 
(FWL) i,i (Fig.2) 

W{S) = jfW{T, -F„ I a e plaquettes of S) = 



^W{r) n W{D). (19) 

plaquettes of S 

where 5* is 2-dimensional surface on the lattice with a given boundary F = dS. 
W involves the product of small Wilson loops over all the plaquettes in S as 
well as a big Wilson loop for F. 

If this model has anything to do with QCD it should posses an oper- 
ator like the Wilson loop which obeys an area law in the confining phase. 
(This property need not be related to existence of quark degrees of freedom). 
Neither the adjoint Wilson Loop nor the Baryon Loop will suffice since they 
obey a perimeter rather than an area law even in the standard Wilson theory 
for finite A^. The FWL (or (|23|), to be exact) seems to be the only reasonable 
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counterpart of the conventional Wilson loop in the framework of the KMM. 
There are several arguments to support this suggestion. 

First of all FWL is a gauge-invariant operator (unhke W{r) which is not 
invariant under local Z[N)). It can thus have a non- vanishing average value. 
Secondly, in the naive continuum limit the FWL reduces to the conventional 
Wilson Loop 

yV{S) r^trPexpi £A^{x)dx^. (20) 

(In particular, it depends only on F = dS.) Indeed, ii U — exp{iaA) then 
the plaquette loop 

W(n) = tr(I + ia'^F - ^a^F^ + ...) = 

= 7V-ia^F2 + ..., (21) 

where F is the field strength. The last term vanishes in the naive continuum 
hmit and thus l^(n) tends to a field independent constant. (Note that we 
have used the fact that tvF = in this derivation). Therefore in this naive 
limit the plaquette loops do not contribute at all. Their only role is to make 
the entire quantity gauge invariant. (Note that along with the ansatz U = 
ex.p{iaA) one should consider a whole class of ansatze U — uj{x,y) exp(mA) 
with an arbitrary oj{x,y) e Z{N) at each link. This would lead to the 
cancellation of all contributions to ^ iy(r) 3>, but not to the average of 
the Filled Wilson Loop <^ W(S') ^.) We shall further justify the use of the 
FWL by arguing later in the paper that the expectation value of the FWL 
obeys the area law, 

<>V(5)> ~exp(-Q;'A(5)), (22) 

where A{S) stands for the area of S (i.e. the number of elementary plaquettes 
in S). 

Before proceeding we must discuss the normalization of the FWL. In the 
Wilson theory every loop l^(r) is usually accompanied by a factor of 
According to this rule the FWL operator should contain an additional factor 
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]\[-MS) jj^ ^]^g definition (0). This would lead to cancellation of all the 
factors of in the naive continuum limit leaving precisely a conventional 
Wilson Loop with its correct normalization factor Notice however that 
we have used a different normalization factor in ([191) - only a single factor 
of j^. The extra factor of makes a very big difference in that it can 

lead to a spurious area law behaviour for the FWL. We shall argue below 
that this normalization is indeed the proper normalization. We shall 
see that the naive continuum limit is unrelated to the true continuum limit 
of the theory and that only this normalization leads to finite averages 
in the true continuum limit. 

An additional complication which is introduced by the FWL and is un- 
avoidable in the KMM is the dependence of its expectation value on the 
surface S rather than just on its boundary F. It is thus nearly certain that 
the correct "physical" operator (i.e. the operator which is to be compared 
with the Wilson loop of continuum QCD) should contain an additional sum 
over all surfaces S with the same boundary F: 

W(F) = f [VS] W{S) (23) 
Js,ds=r 

with some (string-theory-like) measure [T>S]. This is important, for example, 
for the study of correlators of FWL operators (see Section 6). In this paper 
however we restrict ourselves to the problem of the evaluation of -C W(S') 3> 
for a given surface S. (In particular in order to derive the area law in the 
KMM it seems enough to take for S the surface of minimal area with the 
given boundary F.) 

3 Integrals over unitary matrices 

As discussed in the previous section we are interested in evaluating the av- 
erages of operators such as filled Wilson loops. First note that averages of 
any operators which are invariant under the Z{N) gauge transformation (|l^) 
are consistent with the integration over angular variables which is implicit 
in the transformation of ([l|) into (0). Thus the the same operators should be 
averaged in the "Matrix Model" approach (0) as are averaged in the original 
formulation (|l|) of the KMM. To evaluate these averages one must first study 
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integrals over unitary matrices f/ on a single link which are of the form 

= J [dU] expS».. '^''^^l^'^l' ..f/,^,-^f/t'=i'i...f/t'=«'^ (24) 

Due to the Z{N) invariance under U uU, W u'^U'^ the integral 
is non-vanishing only when the number of operators U and under the 
integral differ by an integral multiple of i.e. p — q = [mod N]. (This is 
the case for U G SU{N). li U e U{N) the condition is that p = q). 

The integrals ( p^ are well known for p=q= 0. In fact there is an explicit 
formula 

r det • • e'^^'^^ 

/oo[0,^] = y [rff/]exp(5]0,V',|f/^,n = const^^^^^^. (25) 

The standard references for this result are [^]. There are, however, two more 
contemporary derivations which deserve mentioning. 

One such derivation starts from an interpretation of /qo as an integral 
over the co-adjoint orbit of 

/oo[0,^] = ^/[c^X]e*'-*^ (26) 

where X = U'^U'^ is a generic element of the orbit of \E' under the action 
of G/H, where H is the Cartan subgroup of G and [dX] is the standard 
symplectic measure on this orbit. (Note that the diagonal components of U 
do not act on a diagonal matrix \Ef.) This integral was explicitly evaluated 
in with the use of the Gelfand-Tseytlin parametrization of the orbit. 

Another modern derivation of Iqo is based on application of the Duister- 
maat-Heckmann (DH) theorem The idea of this theorem is that (under 
certain very restrictive conditions) an integral can be substituted by a sum 
of quasi-classical contributions evaluated in the vicinities of all the extrema 
(not only the minima) of the action. Applicability of the DH theorem to /qo 
relies on the fact that /qo is essentially the integral over an orbit which is a 
symplectic and symmetric manifold. 

To see how the DH approach leads to the formula (pSD note that each 
extrema of the action J^ij 4'i'^j\Uij\'^ as a function of U is given by a product 
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of any diagonal matrix Ud and any matrix of permutations P i.e. U = UdXP. 
Diagonal matrices do not contribute to the action, and thus the sum in the 
DH formula goes over all permutations of N natural numbers l...A^. In order 
to evaluate the contribution of any given extremum one needs to find the 
determinant of quadratic fluctuations about this extremum. To do this we 
write U = Pe*^ where if is a Hermitean matrix. The action is then expanded 
to quadratic order in H: 



Then DH theorem implies that all the higher order terms can be ignored. 
The result of the Gaussian integration over H is 



J[dU]explY.<f'^^J\U^J\' 



p -^-^ i,j 



p 



A(0)A(^) A(0)A(^) 



where the sum is over permutations P and the factor (— )^ comes from the 
fact that 

n(V'p» - = i-y Ui^^ - ^.) = i-fH^) (29) 

i<j i<j 



This completes the derivation of the result ( ^31) using the DH theorem. 

Unfortunately the same methods are not directly applicable to evaluation 
of the integrals Ipp with p ^ 0. In fact the generating functional for all non- 
vanishing integrals of this type is given by 

I[A] = J [dU] exp J2 Ai|f/.,f . (30) 
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Unfortunately the DH formula is not directly applicable to this case. If it 
were applicable then we would find that: 



with order 1/A corrections vanishing. (Note that the classical equations of 
motion are J2j Uij{Aji — Ajk)Ujf^ = 0. Any permutation matrix f/ = P is a 
solution to these equations.) Since the DH theorem is not directly applicable 
in this situation the C(^) corrections do not vanish in general. There are 
several interesting exceptions in which case they do vanish. The first case 
is if Aij is of rank 1 so that Aij = (piipj (we then recognize Eq.(p8D in Eq. 
(|3lD). Another interesting case is when N = 2. For the case of SU{2) one 
can check explicitly that Eq. (^) is exact so that 

(Note that for N = 2 the denominator is independent of P). Other excep- 
tional cases involve certain ^/-integrals arising in the study of correlators in 
multi-matrix models like 

''c/$rf^e*'-**tre"*+'^* (33) 

The fact that the DH theorem is not directly applicable in the general case 
( |30[ ) is rather obvious. The action in the integral (|30| ) is not defined on any co- 
adjoint orbit but on the entire group manifold which is not even symplectic. 
Keeping this in mind it is not too difficult to work out the corrected formula 
(see p) though we shall not use it in this paper. In what follows we shall 
make use of the approximate formula ( PTf ) and of the exphcit expression (^) 
for SU{2). 

For our purpose which is the evaluation of the expectation value of the 
FWL we are mostly interested in the correlators Ipp with p = 1: 

If. =< f/.^.[/tfc' >= 1 \d\J\ eS..^^'^^l^-l' [/,^.f/t^'. (34) 
It is useful to use a somewhat more explicit representation of the Haar mea- 
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sure:[] 

[dU] = n dU^.SiY^ U^kWk - (35) 

i,j k 

Making use of the invariance of the measure under the transformations 

Uij^JUij, u^k^ ^U^^^uj^^ (36) 

and 

Ui, ^ Ui.uj^, U^''' uj"''U^^' (37) 

^ This representation of the Haar measure suggests a possible approach to the evaluation 
of the integrals (|30|). The idea is to introduce Lagrange multipliers Xij — Xji for every 
one of the (5- functions in (p5|). Then Gaussian integral over unconstrained [/-variables is 
easily evaluated with the result 



I[A] 



For N — 2 the integral over (PXi2 gives 



/ 



Aii(A21 - A22) + A22(Aii - A12) + {A12A22 - ^11^21) 



n All — All , , A22 — ^211 



All — A12 A22 — ^22 

Consider the first logarithmic term. The integral over A22 is easily evaluated yielding 

1 ^11^21 - A12A22 
I[A\ (X — exp X 

An- A12 An - A12 



f t . An+ A22- A12- A21.. 

/ dAiijexpAii 1 log 

J An — A12 



All - An 



All — A12 

This integral can be evaluated with the help of the formula 

dAe^^ log ^ = - / dXe^\-^ = \e^" - e^^) 

X ~ p c J X — a A — p c 

and is equal to (|3l]). The contribution of the second logarithmic term is is the same as 
the first. 
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(with uj'^ = 1), one can concludes that I^j = unless i = I and k = j, so 
that 

< UijU^''^ >~ Sis';. (38) 
We now introduce the quantities 

a.^ ^^^f^^ , tJ = l...N, (39) 

which are the only non-vanishing quadratic correlators. These clearly satisfy 
the relations 

Q,>0; = 1 = ^C,, forVz; ^] = (40) 

j j 

It thus follows that the Cij can be interpreted as conditional probabilities. 
These quantities will be used in the next section to evaluated the average of 
the FWL. 



4 Loop averages as partition functions of sta- 
tistical systems 

Originally the average -C W(Ti...Tk) ^ for Fi + ... + = can be repre- 
sented as a double-loop diagram ("fat graph"). Because of the property (^) 
all double lines can be exchanged for single lines (Fig. 3) since 

^Iff^ = Su6,,Q,, (41) 

Thus at every site there is only one independent value i of some effective 
"spin". These single lines are oriented because, in general, Cij ^ Cji. Let 
us denote the graph consisting of these single lines by F. For FWL T = S. 
Then 

« W^(Fi...Ffc) » = 



J d^exp[-V{(p)] < 1 > Et{x){U<^.,y>(.rCi^x)jiy)[(p{x),^{y)]} 
Jd(f)exp[-V{<f))] < 1 > 



(42) 
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where the sum is over distributions i{x) of integers i at sites x. The quanti- 
ties on the r.h.s. are aheady very similar to the partition function of some 
statistical system with only nearest neighbor interactions. In the case of the 
FWL operators this system is essentially 2-dimensional (for adjoint loop it 
would be 1-dimensional). It is unfortunately not a very simple theory, since 
the statistical weights at any link depend not only on the values of the "spin" 
at two adjacent sites, but also on the orientation of the link {Cij 7^ Cji) and 
on the field We can thus think of the statistical weights Cij as random 

variables with some non-trivial distribution which is encoded in terms of some 
"hidden variables" (j){x) with gaussian distribution when V is quadratic). Be- 
cause of the presence of these random variables the system in question looks 
more like a "spin glass" than an ordinary spin system. Q 

A drastic simplification arises in the "mean- (or master-) field" approxi- 
mation for the integral over 0(x) in which case (f){x) is frozen and does not 
depend on x, i.e. 0(x) = $. The C^j are thus also fixed and, moreover, 
they are symmetric i.e. Cij = Cij[$, $] = Cji[^, $] = Cji. W{T) is thus the 
partition function of a magnetic-type statistical system: 

« w{t) » ~ E { n c^^^^mv)} (43) 

<x,y>£r 

The fascinating thing about this relationship is that the area law for the 
FWL is just equivalent to the statement that the free energy of the statistical 
system is proportional to its area. This will be true provided that the free 
energy of this 2-dimensional system is nonzero. This seems to be a rather 
obvious property for a system with local (nearest neighbor) interactions. 
There is however one delicate point which prevents us from using this result 
immediately. The problem is that the "Boltzman factors" have the rather 

^ Because of the presence of the factor < 1 > m (^Tj) which is (p dependent, the shape 
of the random distribution of the Cij depends on the actual choice of contour T, i.e. on 
the form of the 2-dimensional sample S = T. This is a kind of a boundary effect for the 
statistical system and although it is is not very important for the derivation of results such 
as the area law it is important for the study of correlators of (filled) Wilson loops. 

Note also that the r.h.s. of Eq.(p^) is a ratio of averages over ^-fields. Formally this 
is different from the situation in conventional spin glass systems where any back reaction 
of the system on the shape of the random "noise" distribution is neglected and thus 
the average of the ratio rather than the ratio of the average is the quantity of interest. 
(In practice the difference is not so drastic, since the distribution of the logarithms of 
Boltzmann weights logCy in ( ^ ) is far from being Gaussian.) 
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peculiar normalization (^Ol). We now discuss the implications of this fact 
using the simplest example of the Ising model which arises for N = 2. 



5 The case of = 2: The Ising model 

When N = 2 the diagonalized matrices and ip can be written as 

^■-f! (44) 




It is now easy to express the Cij explicitly in terms of the variable 7 = 2(f)'ip 
One finds that Q 

_ /o^(l - x)e-'^-dx 1 - 1/27 + e-^V27 

\lxe-^~^^dx 1/27-6-^^(1 + 1/27) 
- - \le--^^dx - ^^^^ 

The case N=2 is mostly distinguished by the fact, that the two diagonal 
element Ca coincide i.e. C\\ = C22 (and, as a result, C12 = C21) for any pair 
{(pjip)- This allows one to represent Cij as 



Cn — C22 



.JIT + e-j/T 

-J IT 



C12 - C21 - gj/T ^ g-J/T (46) 

SO that <^ W{S) ^ is proportional to the partition function Zi^ing of a 
2-dimensional Ising model defined on the surface 5*. In the mean-field ap- 
proximation when (f){x) is constant, J is the same for all sites and this is 

^ These formulas can be easily deduced either by an explicit evaluation of the integrals 
over SU (2) — S'^ or by an application of the DH result ( p2|). The simplest way to proceed 
with an explicit derivation is to use the representation's^) for the measure [dU]. One 
writes U — ao + ia ■ a. The (5-functions in (^5|) then imply that ao and a are real (after a 
[/(I) factor is extracted from the measure), and the measure is simply / d'^aS{a'^ + cP — 1). 
In Eq.(E5|) x = al + al = \Ui2? and 1 - a; = _^ a§ = \Uii\^ = \U22? ■ 
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just a conventional Ising model in a finite area A{S). More generally there 
is some random distribution of the values of J . This distribution will not, 
in general, be Gaussian unless T and I/7 are small. The region of small T 
will be of great interest to us later in this paper. In what follows we restrict 
ourselves to the mean-field approximation (and, without loss of generality we 
set J = 1). 

The average of the FWL can now be evaluated with the help of Eq. (^) 

The two dimensional Ising model is exactly soluble. We can thus substi- 
tute the well known formula for the free energy of the Ising model (see, for 
example, §141 of 0) to find the string tension 

- log < W > 

" = MS) = 



(1 + a;2)2 — 2x(l — a;2)(cosu;i + cosco'2) 
duJiduj2 log (4j 



where 



X = tanh (1/T) 

The three most interesting temperatures in the Ising model are T = 
00, T = Tcrit = 2/ log(A/2 -|- 1) and T = 0. a' is continuously decreasing with 
decreasing temperature T. At T = a'=0 and at T = T^it (Pa' /dT"^ has a 
(logarithmic) singularity. T = 00 corresponds to the value $ = \Jl /'^ = of 

the mean field , while T = corresponds to $ = 00. In terms of the original 
KMM T = cx) is associated with m = 00 (the strong coupling limit), while 
T = Tcrit, where long-range correlations arise in the Ising system, should be 
associated with m = merit (the "weak coupling limit", where continuum- 
like behaviour is supposed to occur). The problem is that at Tcrit there is 
no reason to believe in the validity of the mean field approximation. Our 
main hope is that as increases, the analogue of Tcrit [A^] moves toward T=0 
so that the critical behaviour of the relevant statistical system appears in 
the limit T —>■ 0, (p —>■ 00 and m — > merit at which point the mean field 
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approximation may be trustworthy. It is most difficult of course to prove the 
existence of a second order phase transition for all A^. With this guess in 
mind, let us pay more attention to the point T = in the example of the 
Ising model. 

The crucial feature of the point T=0 for our considerations is that a' — 
as T — > 0. [| It thus looks as if the string tension is vanishing in the continuum 
limit. We should be more careful, however, and consider what happens when 
T is small, but still non-zero. The answer for SU{2) is easily extracted from 



( ^81) from which one finds that 



a' = ^ + o{l/^\e~'^') (49) 

Approximation that is large near the continuum limit is self-consistent. 
The integral over is essentially 



with a' given by (^Uj). Thus when ^crit integral is dominated by 

0^ ~ m? — which is large in this limit. 

We suggest that this simple example of the Ising model where everything 
can be calculated explicitly can serve as a prototype for the situation of 
N = oo. We speculate that the main features of this example survive in the 
general situation. There are, however, several crucial differences: 

(a) The phase transition occurs at a lower temperature Tcrit which ap- 
proaches T = for large N. 

(b) The single eigenvalue gets replaced by a distribution of eigenvalues 
with a density p(0). 

(c) The exp(— 0^)-type corrections can become relevant for smooth func- 
tions p(0). 

^ Note, that this result depends crucially on the normalization of the FWL operator. 
If there were an extra factor of N~^^'^°- = 2~^^^'^ in (p^, we would get a' log2 as 
T ^ 0. This would exclude any dynamics from the derivation of area law which would 
be valid in the continuum limit without any need for scaling arguments. This is already 
an argument in favor of our normalization choice. We shall return to the normalization 
problem in Section 6. 
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If this suggestion is correct then we would expect Eq. (^) to be replaced 

by 



/ / p{<l))p{<l)')d(j)d(f)' 



(note that 20 in (^9]) is just the difference 0i — 02 of two eigenvalues of 
It thus follows that oi is not necessarily vanishing in the continuum limit. 
Critical behaviour would arise at large but finite values of the master field 



6 Puzzles of the continuum limit of KMM 

Before we turn to the general analysis of the S\J (N) case, let us examine in 
more detail the two limiting cases of strong and weak coupling. 

The strong coupling limit of the model corresponds to m — >■ oo. The 
exp — m^0^ factor in the integrand then damps all fluctuations of the fields, 
which are thus restricted to vanish. It follows that the action in the integral 
( |34D over U vanishes and what remains in the evaluation of expectation values 
are just integrals with the conventional Haar measure, of the form 

<i> '^=° - ndu] - iv^"^^-^ (^^^ 

so that 

QA<P = ^ = 0] = ^ (53) 
is equal for all i and j. (One can similarly derive expressions such as 



^ ^ ^ \ij>=0 -^x^il^jk^mq^np ' "iq"jp"ml"nk) 



J\j(^j\j2 '^{^jk^lm^np^qi + ^ jp^ qm^nk^li) ■ (54) 
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In this strong coupling approximation any average of a gauge invariant prod- 
uct of Wilson loops (2.4) is simply equal to 

< W{Ti...Tk) > \strong coupling = j^^Hnks ^ J^#loops ^ jyfc-l (^^) 

provided Fi + ... +Tk = 0. For the FWL operator A; — 1 = 7^ (of plaquettes 
in S) = ^(5'), and we obtain: 

< W(Ti...Tk) > Istrong coupling = j^^rea (^^) 

The area law is thus trivially true in the strong coupling limit. Of course, the 
really interesting question is whether it is also true in the "weak coupling" 
(continuum) limit. 

Naively the "weak coupling limit" corresponds to m ^ 0, where integrals 
over (j) are dominated by large values of the (/)-fields (i.e. all the eigenvalues 
(pi are large). In fact integrals over diverge already for a finite value of m = 
'TT'crit (quasi-classically m'^^-it = D), and the relevant limit is m ^ ?^crit + 0. 
The fact that integrals over diverge still allows interesting quantities such 
as ^ — , which are ratios of divergent integrals, to be well defined. In fact, 
when all (pi's and ipi^^ are large (and positive) the integral 

|[dC/]exp^0,^,|f/,/ (57) 

is clearly dominated by diagonal matrices 

U,,=e''^6,,, |f/^,f = 1 (58) 

and thus 

< UijUli > </,,V)>l jr c c ^ 0,-!/'>l . /_„x 

— — OiidjkOij, Cij —>■ dij (59) 

In this approximation it is clear that all the spins in the statistical system 
take on the same value at each site. (This is clearly the ordered phase which 
occurs in the spin system at zero temperature.) There is just an overall 
degeneracy factor of N, which cancels the only N in the normalization in 

< iy(Fi...F,.) > = 1 X 5(Fi + ... + Tk) (60) 
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In particular in this approximation the FWL is given by 



« W{S) » |„_, = 1 (61) 

and, as explained at the end of the previous section, the vicinity of the point 
m = merit should be investigated in order to explain whether and why a' 
may be non-vanishing in the continuum limit. 

This will be our purpose in the next section. In this section we discuss 
what happens in the "naive continuum limit". As we already mentioned 
in Section 2, evaluation of VV(S') in the naive continuum limit would give a 
factor A^'^"^ instead of 1 (there is fc — 1 instead of k in the exponent is because 
of normalization factor in (|T9|)). This estimate was however based on 
the imphcit suggestion that 

^''^ \ncl = ^ijSkl (62) 

which is obviously inconsistent with the basic property (0) of the correlators 
of U (which can of course be traced back to the Z{N) gauge symmetry). This 
argument is already sufficient to demonstrate the failure of naive continuum 
limit. We shall now present a somewhat more constructive argument against 
the naive continuum limit. 

The discrepancy between the naive continuum limit (which we claim is 
incorrect for the KMM) and the limit of m — merit for the eigenvalue model 
(0) (which we claim is the correct prescription for the KMM) appears already 
in a very simple situation, which can be examined exactly. Indeed consider a 
lattice which is just a closed 1-dimensional chain of L links. All the techniques 
and results which we derived for a generic D-dimensional lattices is applicable 
to this simplified model. The set of (gauge-invariant) observables is very 
small. In fact only adjoint Wilson loops and "Baryon" loops along the closed 
chain and their "powers" (the same contour passed several times) are nonzero. 
Such adjoint Wilson loops are given by 

0, = ^\W{kT)\'. (63) 

The predictions of our above analysis are that in the strong coupling limit 

JY# of vertices \ 
< Ofe > |m=oo = ^ ^# of links = NN^L = i^"^) 
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while in the weak couphng hmit (near merit) 

< CJfc > U=mcrit ^ ^ ■ ^ ^ ■'- in the hmit m merit (65) 
On the other hand the "naive continuum hmit" prescription ( |B2[ ) would give 

< Ofc > \nci = J^M? = (66) 

which differs from (|65D. 

We shall now solve the original KMM (|1]) exactly for this particular lattice 
and prove that while ( |64|) and (|65D are correct predictions the "prediction" 
( |66[ ) is incorrect. We shall also find a reason for the failure of ( |66| ) which can 
serve as a general explanation of why the naive continuum limit is wrong 
in the KMM in higher dimensions. Recall that the argument in support of 
(|66D was that the action in (H) diverges as m — merit and it has a single 
maximum as a function of U which occurs when U = I. This is absolutely 
true but what is unusual is that in the KMM is that the peak of the action at 
U = I is very moderate. In fact the action grows only logarithmically as U 
approaches /. This makes it necessary to take the measure of integration [dU] 
into consideration. The vicinity of U = I where action is big has, of course, 
a small volume and it appears that this is enough to compensate completely 
for the peak of the action. In the particular example of the one dimensional 
chain we shall see that when m = merit the action and the measure cancel 
each other completely over the entire range of integration over U. Once 
observed, this phenomenon is clearly a very general feature of the KMM. 
The integrand of the partition function is very smooth as the continuum 
limit is approached so that even quasi-classically it has is no dominant peak. 
The continuum limit is thus far from naive. 

Let us now turn to our calculation for the one dimensional chain. The 
quantity to compute is 

< Ofe > = 



1 



,x 



+ije 



-m^tr#J 



t 

x,x-\-l 



itr[n u. 



lfc|2 



x,x+l 



(67) 



x=l 



x=l 
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with 



n d^x[dU, 



rt 



X = l 



{x=L + 1 is identified with x=l). This one dimensional case is particularly 
simple since it is possible to perform a gauge transformation ([T3|) which 
eliminates (i.e. makes equal to /) all the matrices Ux,x+i except for a single 
one, say U = Ul,i- In other words 



[dU] n / d^.e 



x=l 



L-1 
x=l 



tr$:j,*^+igtr$z,(7*iC/t 



(68) 



and <^ ^ is simply the normalized average of jtrf/'^p with this generating 
functional. 

The next step is to integrate first over $2, ^l-i and then over $1. All 
these integrals are Gaussian and are thus easy to perform. The result can be 
expressed in terms of the determinant ni of the / x / matrix 



/ 2m2 -1 

-1 2m2 



Ml 










V 



-1 2m^ 



\ 



















-1 2m^ J 



(69) 



The determinant /i; = detMi satisfies the recursion relation /i^+i = 2m? — 
with "initial conditions" /zq = 1, A'-i = 2m^ (in this respect they are 
somewhat similar to Fibonacchi numbers fi+i=fi + fi-i, /o=/i=l). Since 
the matrix Aii is simply a Laplacian the fii can also be represented as 



2'Y[{m' 

k=l 



COS ■ 



/TT' 



(l-g)(l-z/g)---(l 



,21+1 



q) l-q 



21+2 



q\l-q^) 



(70) 



where 



g = g± = ± V m'^ — 1 



(71) 
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and V = vi = e 2'+2 . 

We now do the integral over $^ in (p3): 



■' x=2 •' 

$1 and $L can be considered as "source terms" in the Gaussian integral over 
$2, with the action 

1 

- il^,{ML-2).y% = (73) 

^ x,y=2 

exp ^[$?(A<^i2)22 + 2<^,<t>LiMi'_,h,L-i + <I>1{M,\)l^,,l-i] 
The relevant matrix elements are 

iMl'_,h2 = (-MZ-2)l-1,L-1 = ^ 

and iML'_,)2,L-i = (74) 

/iL-2 



Since - /iL_3/2/iL_2 = /iL_i/2/iL_2 we get for (^: 
/iZ52/^ / rf<|.irf$^exph^tr(<|.? + $i) + 



^ tr$i$L + tr$L[/<l>i[/^]. (75) 



I^L-2 



The remaining integration over $i is done by rescaling — $ = ^l/v^A^l-i 
The result is: 

Z ^ATlJ [dU] J e-^'l''^"" + (76) 

with A/l = 1 and "renormalized" (generically lattice and contour-dependent) 
mass parameter 

^2 _ - /^f-2 - 1 _ /^2L-2 - 1 ^ 1 + 
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In terms of the new parameter qi = Ml — ^ — 1 (i.e. M| = {l + qD /2qL), 
Eq.([77|) turns into a simple relation: 

QL = (g±)^^ (78) 

where q± is defined in ([H]). 

The integral for Z and for <^ (9^. ^ is easily evaluated by performing a 
gauge rotation on U and $ so as to make U diagonal {U = diag(e*^S e*^^)). 
It follows that 

N 

i=i 

while the Haar measure is 

N N N 

[du] ~ n ^'(e'") = n n(e''' - ^''^f ~ 

j=l j=l i<j 

3=1 i<j ^ 

For diagonal \J the action in ( [76D becomes 

TV 



(Ml - 1) 5: + 2 ^ (M2 - cos(^, - e,)) (80) 



The integral over $ is now easily evaluated. Using ([79|) one finds 



1 r"" sm 



2 



2"''^ / n n ^ n ^ _ m 
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The critical value of M is clearly Merit = 1-0 At this critical point the 
integrand in (^) becomes equal to unity. The Jacobian factor in the measure, 
YliKj sin^ -^^2^ completely cancels the contribution from the action, ni<j(l ~ 
cos(^j — 9j)). Note that the contribution of the action is in in denominator 
and is indeed singular as M — > Merit but, as stated above, the singularity is 
very moderate. We may now evaluate the average of |tr?7^p at M = M^it. 
Since the measure of integration is just ]\^ dOj we find 

« O, » U.„„,. - j^^^ - 1 (82) 



in accordance with (^) (and contrary to (|66D). This completes the proof. 
The conclusions which we deduce from this result have been discussed pre- 
viously. 

In the remainder of this section we use the opportunity provided by this 
calculation to make another important observation. Note, that there is a 
somewhat different way to evaluate the integral ([76|) . It can be represented 
as the determinant of an x A^^ matrix, 

(76) ~ det-^/^{M2/(g)J--(f/®[/t + [/t®f/)} ~ 

2 2iW"^ 
M /^exp- > - — -7tr( )' 



^Eq. ([ttI) implies that M^j.;j = 1 (gcrit — 1) simply corresponds to^j,;j = 1, which is the 
same for all values of L. (This is easy to check, since for m? = 1 the formula for m is very 
simple: /i; = l + \ This value of m^j,jt is also in accordance with the simple quasi-classical 
estimate tti^^jj = D (see footnote 3 above). 

Eq.(|8l|) was also deduced in and in a very recent paper [jll], where also the result 
of integration over all ^-variables was derived: 

' (HI) - T, ^(T^^ n —y 9l = fe)^"^- 
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If we now substitute U in the diagonal form this expression becomes 
M-^V2 exp i V -i- V ±^ y = 



^-''-^^jtl^ ^ tt}C-^)^} ^ (84) 



1 1 1 



i.e. we reproduce the contribution of the action to 



The exponent at the r.h.s. of (|8^) can now be rewritten as 

1 °° 1 ' /' 

2 S w £ ^w^'''^'"' ^^^^ 

In this expression one can easily recognize the D = 1 version of (|^). 

We now review our conclusions from this section. Our main claim is 
that the action grows only logarithmically in the vicinity of the critical point 
m = merit- This fact is of course related to the occurrence of the factor l/l 
in the sum over Wilson loops in Eqs.(^,^). However, the paths which are 
summed over in these formulae are allowed to have steps which "double back" 
on themselves (Fig. 4). Sums over such inverse steps introduce comphcated 
combinatorial factors such as the binomial coefficients in (pS]). The above 
calculation serves as an illustration that these complicated corrections do not 
necessarily destroy the logarithmic behaviour but rather renormalize the bare 
mass parameter. 



7 The area law in the case of arbitrary N 

So far we have examined two exactly solvable examples. We have solved the 
KMM for the group SU{2) and for 1-dimensional lattices. We return now 
to the general case where N is arbitrary and the lattice is D-dimensional for 
which an exact solution is not yet available. We now suggest a qualitative 
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picture which we expect to be vahd in the generic situation. This would 
include the following ingredients. 

(a) Gauge invariant observables in the KMM are associated with either 
one or two dimensional sublattices F of original lattice. The corresponding 
average is given by the partition function of a Id or of a 2d statistical system 
on r (which is defined by the Boltzmann weights Cij[(j), ip] with randomly dis- 
tributed 0, ip) . FWL operators are associated with two dimensional lattices 
r and thus with two dimensional statistical systems. 

(b) The area law for the average of the FWL is simply a reflection of the 
fact that the free energy of a 2d statistical system is proportional to its size 
(area). 

(c) In the strong coupling limit m — > oo the string tension a' = log A^. 

(d) Long range correlations occur in 2d statistical system, if they have a 
second order phase transition, at a certain temperature T[N;m\ = Tcrit[N] 
corresponding to a certain value of m = merit [A^]- It seems that a necessary 
condition for the KMM to have a continuum limit is the presence of long 
range correlations in the 2d subsystem and thus the phase transition at Tcrit 
must be second order. 

(e) The partition function of the 2d statistical system can easily be 
evaluated in the vicinity of an a-priori different point r[A^; m] = when 
m = merit [-/V]- (The partition function of the KMM diverges at this point. 
This is why the subscript "crit" seems reasonable.) When T is precisely equal 
to zero the string tension a' = 0. However, in the vicinity of the point T = 
it obeys some scaling law such as a' ~ l/a?N"' . It can thus be finite in the 
appropriate limit a — 0, — oo. The region of small temperatures is easy 
to analyze since fluctuations of the Boltzmann weights are suppressed i.e. 
the 0-fields are "frozen" and the statistical system is in its "ordered" phase. 

(f) It may be hard to find the critical temperature Tcrit [iV] for the actual 
system (|2]). If, however, (^) is evaluated in the (unjustifiable) approxima- 
tion of frozen 0-fields then Tcrit[iV] > 0, and, moreover, the phase transition 
is generically first order. However, in the limit N ^ oo both Tcrit[iV] — 
and the latent heat A — >■ 0. The reason for this is that the number of degrees 
of freedom and thus the fluctuations are increasing as — oo. This increase 
in the fluctuations should also occur for finite A^ < oo if one goes beyond 
the frozen-0 approximation. Thus in the large A^ limit the two interesting 
values of T namely T = T^rit from d) and T = from e) seem to approach 
each other. Thus the continuum limit of the KMM (if it exists at all) can 
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be examined in terms of the phase transitions of 2d statistical systems. We 
thus expect the continuum hmit to occur at small "temperatures" . 

(g) Our description of the continuum limit is somewhat more intricate 
than the usual description. This is a peculiarity of the KMM where the "naive 
continuum limit" which is used in the Wilson theory (/[c?t/] exp(^ J2 W{0)-\- 

c.c.)) and in the "adjoint" theory {J[dU] exp(^ J2 {without a sum 

over non-elementary contours!)) does not exist. Although the most natural 
assumption after the failure of the "naive continuum limit" would be that 
there is no continuum limit at all, the above scenario implies that it still can 
exist but that the approach to this limit is more sophisticated and interesting 
in that it is described in terms of 2c? systems! 

The remainder of this section is devoted to various comments on this 
scenario. Most of them concern notes d)-f ) above. The note g) was already 
discussed in some detail (though not exhaustively) in Section 6. Notes a)-c) 
were discussed in Sections 2 and 4 and illustrated by the example of the Ising 
model in Section 5. Our goal now is to try to generalize some results of that 
section to iV > 2. 

We begin by deriving the literal analogue of Eq. (^9|) for the average of the 
FWL for the case of generic N. We shall subsequently discuss the validity of 
the result. For this derivation we simply need to find the first order correction 
to the free energy of our statistical system at low temperature in the master- 
field approximation for assuming that the system is in its ordered phase. 
Since at T = Cij = 6ij we introduce a new variable ^ij where Cij = Sij 
In the ordered phase the "spins" at all the sites of the sublattice T = S have 
the same value i. Taking into account the elementary fluctuations (i.e. the 
flip of a "spin" i j at one point of the lattice) as well as the deviation of 
Cii from unity, we get: 



Due to the normalization condition (|y) = —J^j^ti^ij, while the second 



a 



log < WjS) > 
a^A{S) 




i^ii + y^AiF) + higher order corrections (86) 
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term in (^) can be neglected since 2D > 1. Thus 

a'^^Efc. (87) 

In order to estimate in the limit of large (which is the same as "low 
temperature" and m —> merit), one can make use of the formula (|3l|) (C^(^) 
corrections being unessential in this limit). 

The leading contribution to (|3T|) comes from the identity permutation 
P = I. The corresponding contribution to is 



d , e^k ^fefe / /I 



dAij ^ l\k<iiAkk + All - Aki - Aik) \ \A 
Since we are only interested in the case i ^ j, this gives 

0(i 



An + Ajj — Aij — Aji \ \A 



2 



1 + 0{h] (89) 



In the approximation where all the 0's are large and different, the O (1/0)- 
corrections can be neglected. The same is true about the contributions of 
other permutations P I, which are suppressed exponentially: the contri- 
bution of Pki, which permutes k and /, is O (^e~^'f'''~'^''^^^. Therefore the literal 
analogue of (^9|) for generic N is 



However as ^ oo the most interesting phase, associated with the 
continuum limit of the KMM, seems to occur when eigenvalues (pi, though 
large (0i ^ 1), are not all very different i.e. ^1. In other 

words the eigenvalues should be smoothly distributed with some density p(0) 
(/ p{'t>)d'(p = N), so that Y,iF{(j)i) ~ / p{(j))d(j)F{(f)). If expressed in terms of 
p{(l)), Eq.dig) looks like 

2 ,p(^J^(^ 
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For smooth distributions p(0) the corrections, neglected in (p8D, can be- 
come essential. First of all there can be corrections which contain negative 
powers of the differences (0j — (pj) which are enhanced in integrals like (PT|). 



Indeed the (9(l/A)-corrections in (^) contain terms with negative powers 

(Ap(i)+^jPO)-Ap(i)-^jP«)~'' ^ with k = l- ■■Tank{SU{N)). 

If it were not the k = 1 term (as supposed in (|8^)) but the k = N — 1 term 
which dominated then we would get an extra factor of (A^ — 1) on the r.h.s. 
of (pOl). It may however happen that the contributions of all k are equally 
important in which case the exact counterpart of (|3TD should be used (see 
[|]). Secondly when {(pk — fpiY is not usually large the exponential corrections, 
(arising from non-trivial permutations P), are also important. Expressed in 
terms of the density p(0) they have the form 

3("^-^')'G(0, (j)')p{(f))p{(f)')d(j)d(j)' (92) 

where G is some rational function of and 0'. 

We shall not go into more detail about the evaluation of these correc- 
tions because the very approximation (the low temperature expansion for 
the ordered phase) may be unreliable in the situation of interest. As we have 
stated previously, in order for a non-trivial continuum limit to occur there 
are several different phenomena which must occur simultaneously. 

(a) The auxiliary 2-dimensional statistical system must have long range 
correlations. This is most probably realized at the critical point Tcrit which 
should be (nearly) second order. 

(b) The master-field approximation for the 0-fields must be applicable. 
This is a reasonable approximation when N oo and it leads to the occur- 
rence of magnetic-type statistical systems. 

(c) The "temperature" T which characterizes the behaviour of the 2- 
dimensional system should be close to zero. At T = the string tension van- 
ishes, a' = 0. This is a normalization condition for the "Boltzman weights" 
Cij related to (^Ol). For the statistical system it corresponds to subtracting 
the ground state energy and is independent of other details of the system. 
On the other hand the area law which occurs at small but non-vanishing 
temperature depends on the the actual dynamics of the statistical system, 
i.e. on whether there is a scaling limit at T ~ 0. 

These requirements are mutually consistent at large A^, when TcritfA^] is 
assumed to approach zero. This in turn implies that T = is in fact a 
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critical point and thus one cannot rely on the low temperature expansion in 
the vicinity of this point (as we have done for the Ising model). One should 
go beyond this approximation and apply some formalism adequate for the 
study of phase transitions. This should be, however, a solvable problem, since 
it concerns a phase transition in a 2-dimensions (though the KMM itself is 
formulated in D dimensions)!. 

It deserves mentioning that in contrast with the situation in the Wilson 
theory nothing special happens to the quantities of interest when we go 
through the (2-dimensional!) critical point. As we saw in the examples 
for finite N the area law for the average of a FWL is equally valid both 
above and below the phase transition: the averages in the KMM are related 
not to correlators of the 2d system, but to extensive quantities (like density 
of free energy). 

It remains to comment on the very notion of "temperature" for these two 
dimensional systems. As is already clear from the example of the Ising model 
in Section 5, this "temperature" is nothing but some function of the "frozen" 
</)-fields. In the case of iV = 2 it is simply related to the Boltzmann weights 
Cij by Eg . (^61) . Denominators in this formula account for the normalization 
condition (^) for the Boltzmann weights, which is peculiar for the KMM 
(and in fact guarantees that a' = for T = 0). This condition cannot 
however be taken into account in the same simple way when N ^ 2 since 
in the latter case the "bond energies" Jij are T-dependent. One can deal 
with this by assuming that the bond energies also depend on some other 
external parameters (like magnetic fields) and the KMM is associated with 
some particular line in the space of all of these parameters, labeled by T. 
This remark is intended to alleviate confusion resulting from the counter- 
intuitive properties of the system which arise from the normalization of the 
Boltzman weights. The constraints (^) guarantee that they are conditional 
probabilities which is unusual for Gibbs-like distributions. 

Thus we see, that much remains to be done in order to study the contin- 
uum limit of the KMM, both for finite (where it may also exist!) and for 
infinite A^. In this paper we have tried to show that this study is very much 
in line with modern mathematical physics and it should deal with familiar 
2-dimensional systems, their critical behaviour (i.e. conformal models and 
their perturbations) and string-like models. 
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8 Conclusion 



In our preliminary investigation of gauge invariant observables in the KMM 
(|l]) we find tlie model to be very interesting on its own (irrespective of whether 
it is related to 4-dimensional gluodynamics). It is related in a non-trivial 
fashion to various statistical systems and a deeper understanding of these 
connections seems desirable. Moreover the study of the KMM for finite 
values of N seems to be already giving us non-trivial information about the 
probable large behaviour - a fact which is rather common in modern 
approaches to matrix models. 

As to the relation of the KMM to gluodynamics we emphasize primarily 
the drastic difference between the behaviour of the KMM as the continuum 
limit is approached and that of the conventional Wilson lattice QCD. The 
KMM approaches this limit in a much smoother way. This may be an ad- 
vantage of the model, at least from the viewpoint of exact solvability. It 
however implies that the transition to the continuum limit is far less naive 
than is usually believed. (We illustrated this point with an example of the 
normalization of observable operators.) This phenomenon is related to the 
relevance of the angular variables (rather than only the eigenvalues) of the 
$ field in the KMM and it can thus be of general importance for all c > 1 
models. 

Another difference between the KMM and the Wilson theory - the oc- 
currence of extra Z{N) gauge symmetry - may somehow be related to the 
previous difference. In fact whenever this Z{N) symmetry is unbroken it 
comes very close to guaranteeing the existence of a confinement phase all by 
itself since all operators with any color properties, even locally, are simply 
not gauge invariant. The operators which can be considered look like nets 
of double lines imitating planar-diagram structures. Despite this essentially 
confinement-like description, the area law in the continuum limit seems to 
arise in a somewhat non-trivial manner which accounts for the actual dy- 
namics of the theory. 

We argued that the dynamical behaviour of the KMM can be studied in 
terms of 2-di'mensional systems. This fact reflects the already mentioned 
difference between the KMM (and, in fact, the "intermediate" model (pi])) 
and the Wilson theory which makes KMM especially attractive, namely that 
the string-like description is essentially built into the theory. We have already 
seen this in several places. Elementary gauge-invariant observables like the 
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FWL depend on a surface rather than just on a contour. Moreover, in 
reasonable approximations (hke that of the master-field approximation for 
0) their averages can be evaluated within the framework of auxiliary 2 — 
dimensional models. 

The operators related to the FWL which are presumably the physically 
relevant ones were suggested in ( ^3]) to be of the form >V(r) = J[DS]yV{S) 
and thus to include sums over all surfaces with a given boundary. Note, 
that the averaging procedure, ^ prescribed by (|l|) does not involve 

any summation of this kind. It comes entirely from the definition of the ob- 
servables. Conversely the "stringy" measure [dS] need not involve the usual 
g-Area fg^^^Qj- gjnce it appears from the area law behaviour of -C W(S') 

If we now consider the correlator of two Wilson loops the quantity of 
interest is 

G= <w(ri;r2)> - <w(ri)> <w(r2)> = 

/ [DS] < >V(5) > - (93) 

JdS=ri+r2 

[ [VSi] I [VS2] < W(5i) > < W{S2) > . 

The first term on the r.h.s. contains sums over connected and disconnected 
surfaces (Fig. 5). and 

G = < >V(ri; T2) >connected + 

< w(ri)>v(r2) > - < w(ri) > < w(r2) > (94) 

The last two terms on the r.h.s. are equal to 
J [DSi][VS2] (< >V(^i)>V(^2) > - < >V(^i) > < >V(^2) >) . (95) 

The expression in the brackets obviously vanishes in the master-field approx- 
imation for (f) unless the surfaces 5*1 and 5*2 are intersecting. However as soon 
as D > A two 2-dimensional surfaces in general locations do not intersect at 
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all. For D = 4 they intersect only at isolated points. It follows that discon- 
nected surfaces do not contribute to the correlators (perhaps, up to contact 
terms), which are thus represented as sums over connected 2-dimensional 
surfaces with a fixed boundary distributed according to the area law. (At 
large distances the minimal-area surfaces will dominate.) 

Another example of the same type arises when the KMM theory is con- 
sidered at finite temperature (not to be confused with the "temperature" 
in our analysis of auxiliary 2-dimensional systems). What is peculiar to this 
situation in Yang- Mills theory is the occurrence of uncontractable closed con- 
tours C which wrap around the periodic imaginary time r with a period /3. 
This contour has an associated Polyakov loop (PL) operator: 

L{C) = WiC) tiPexpifAodr (96) 

<^ 

Usually this operator is the relevant order parameter for the confinement- 
deconfinement phase transition. It is however not suitable for KMM since 
it is not invariant under the Z{N) gauge symmetry so that -C L ^= for 
all temperatures. This is, of course, contrary to the usual case where in the 
high-temperature phase <^ L 3>7^ 0, signaling the breakdown of a global Z^- 
symmetry. Our standard trick of "filling" the Polyakov loop is not applicable 
since the contour in uncontractable. What can be considered instead in the 
KMM is the counterpart of a correlator of two Polyakov loops (Fig.6): 

<i:(c,c")> =L{c)L{c') n ^(°) (97) 

□e5(c,c") 

One can consider the behaviour of this "Filled Polyakov Loop" (FPL) when 
the distance R between C and C is large. In analogy with Polyakov's result 
for the ordinary QCD one can suggest that there are two possible regimes: 

«/:(C,C')»^e-"('^)^; /3>/5„it 
< C{C, C) const] (3 < pcrit (98) 

and there is a confinement deconfinement transition at some critical inverse 
temperature /3crit, such that m(/3crit) = 0. 

^ What is amusing is that precisely the same configuration was suggested in to 
describe the vortex-antivortex configuration in the matrix model for the d = 1 string. One 
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Thus, if one accepts that observables in KMM can be represented as 
sums over fluctuating 2d surfaces, then their contractable boundaries (Wilson 
Loops) play the role of the physical states (gluonia), while uncontractable 
boundaries (Polyakov Lines) describe (world-sheet) vortices. We see, that 
in contrast to the Wilson theory the string-like description is somehow built 
into the KMM. It just cannot be treated in any other terms. This may be the 
main reason why this model should be studied in more detail and why one 
may ultimately believe in its relevance for the theoretical analysis of realistic 
QCD. 
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can consider the boundary Polyakov lines L{C) and L^(C") as a vortex and antivortex on 
the world sheet given by the net nnes(c C) This coincidence is not accidental and 

is connected to the fact that the confinement-deconfinement (or Hagcdorn ) transition is 
equivalent to the Berezinski-Kosterlitz-Thouless (BKT) transition on the world-sheet 
which involves condensation of the vortex-antivortex pairs. The critical phase transition 
occurs when the vortex-antivortex configurations become unsuppressed, i.e. precisely when 
^ £(C,C') ^ 7^ 0. Analogy between the confinement-deconfineinent transition in large 
N QCD and the BKT transition on the world sheet was also considered in p^ . 
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Figure Captions: 



Fig. 1 Adjoint Wilson Loop for a rectangular contour T. 

Fig. 2 Filled Wilson Loop for a rectangular surface S. 

Fig. 3 Example of the substitution of the "fat graph" Fi + r2 + 
...Fg = by an ordinary graph F according to the identity (|41|) . 



Fig. 4 Examples of contours with "inverse steps" which are 
included in the sums @ and ( ^) . Only contours of the type (b) 
exist in the case when the lattice is a one dimensional chain. 

Fig. 5 The identity < W(Fi;F2) >=< W(Fi,F2) >connectcd + < 
W(Fi)W(F2) >. 

Fig. 6 The surface S arising in the definition of the filled Polyakov 
Line C{C,C'). 
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